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In Quantum Chromodynamics (QCD) the eigenmodes of the Dirac operator with small absolute 
eigenvalues have a close relationship to the dynamical breaking of the chiral symmetry. In a sim- 
ulation with two dynamical quarks, we study the behavior of meson propagators when removing 
increasingly more of those modes in the valence sector, thus partially removing effects of chiral 
symmetry breaking. We find that some of the symmetry aspects are restored (e.g., the masses of p 
and fli approach each other) while confining properties persist. 

PACS numbers: ll.15.Ha, 12.38.Gc 
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I. MOTIVATION AND INTRODUCTION 

Dynamical chiral symmetry breaking in QCD is asso- 
ciated with the low lying spectral modes of the Dirac 
operator D jlj. They affect the path integral weight of 
the gauge configurations through the determinant of D. 
As indicated by the Atiyah-Singer index theorem [2| , the 
exact zero modes are related to topological excitations, 
the instantons. For Dirac operators violating chiral sym- 
metry these are real eigenmodes. The nearby non-real 
modes are also thought to be related to composed struc- 
ture of, e.g., overlapping instantons 3]|j 

In a series of papers ^^ it was emphasized that low 
modes saturate the pseudoscalar and axial vector cor- 
relators at large distances and do not affect the part 
where high-lying states appear. In [0, 01 low mode sat- 
uration and also effects of low mode removal for mesons 
were studied for quenched configurations with the over- 
lap Dirac operator [8|, Q • 

Subsequently low modes were utilized to improve the 
convergence of the determination of hadron propagators 
[a, 0, [ifl \lM (s6s ^Iso the recent study [IZ, 13] compar- 
ing the efficiency when using the low modes of the Dirac 
operator or the hermitian Dirac operator, where strong 
dependence on the parity of the hadron states was pre- 
sented) . 

Associating the low mode sector with the nonpertur- 
bative chiral symmetry breaking and the condensate [l|, 
a complementary question is how important it is for con- 
finement and mass generation of hadrons. Here we study 
what happens, if one removes up to 512 low lying modes 
from the valence quark sector. We compute propagators 
of the pion and other mesons and determine the effect of 
this removal on the mass spectrum. This way we want 
to shed light on the role of the condensate related to 



the spectral part of the Dirac operator in confinement 
and chiral symmetry breaking. Our analysis is done for 
configurations generated for two light, mass degenerate 
dynamical quark flavors. The removal of the low lying 
modes is effective only in the valence quarks sector. How- 
ever, as will be seen, this already has significant impact 
on the meson mass spectrum. 

In l4|l5| it has been conjectured that chiral symmetry 
is "effectively restored" for highly excited hadrons, in 
the sense that valence quarks become less affected by 
the quark condensate. This situation is similar to ours, 
where we artificially suppress the condensate as seen by 
the valence quarks. In the context of effective restoration 
such an approach has been discussed already in [y, Il6| . 



II. REDUCED DIRAC OPERATOR 

Lattice Wilson Dirac operators and approximate 
Ginsparg- Wilson Dirac operators are 75-hermitian, 
IsDjs = D\ but non-normal, thus their spectral rep- 
resentation has real and complex eigenvalues and the left 
and right eigenvectors are bi-orthogonal, i.e. {Li\Rj) = 
6ij . The so-called hermitian Dirac operator D5 = 75 D 
has real eigenvalues fii and the eigenvectors are orthogo- 
nal. 

We want to construct meson correlators from valence 
quark propagators which exclude the lowest part of the 
Dirac spectrum. There are two alternative definitions of 
reduction: based on eigenmodes of D or based on eigen- 
modes of the hermitian Dirac operator. We introduce the 
reduced quark propagator via the spectral representation 
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'^Even when studying the low lying modes in quenched gauge 
ensembles one observes non-vanishing density and also the Gell- 
Mann-Oakes-Renner relation works down to small values of the 
valence quark mass until quenched chiral logs destroy the leading 
chiral symmetry breaking behavior. 



Another alternative works with the bi-orthogonal eigen- 
system of D. The two types of truncation are not equiv- 
alent. We first tested the convergence of the low mode 
approximation and, as has been observed in [13], find 
a clearly slower convergence rate for the standard non- 
hermitian as compared to the hermitian Dirac operator. 



In our study we therefore concentrate on our results from 
truncating the hermitian Dirac operator. 



III. CHIRAL SYMMETRY AND ITS 
BREAKING 

The nonvanishing quark masses of the two Hghtest 
quark flavors are relatively small in comparison to the 
typical QCD scale. Neglecting the masses of the u and d 
quarks the QCD Lagrangian is invariant under the sym- 
metry group 



SU(2)i X SVi2)R X U(l)y X U(l),4 



(2) 



The chiral symmetry SU(2)l x SU(2)i^' consists of in- 
dependent transformations in the isospin space for the 
left- and right-handed quark fields and can be represented 
equivalently by independent isospin and axial rotations 
for the combined quark fields. 

The isospin axial transformation mixes states with op- 
posite parity but the same spin. Depending on quan- 
tum numbers the chiral partners can have the same or 
different isospin. The non-degenerate masses of parity 
partners indicate the dynamical (spontaneous) breaking 
of this chiral symmetry with the order parameter (V'V')j 
the chiral condensate. Spontaneous breaking of the chi- 
ral symmetry leads to the appearance of the pseudoscalar 
Goldstone bosons, the pions. 

The flavor singlet axial transformation symmetry 
\J{1)a is broken explicitly due to the non-invariance 
of the fermion integration measure, the so-called axial 
anomaly. It is not a symmetry of the quantized QCD. 
Consequently no isosinglet Goldstone boson exists within 
the two-flavor QCD and the 77 meson(s) are heavier then 
the pion, attributed to the anomaly. In addition to the 
anomaly also the chiral condensate breaks this symmetry. 

Both symmetry breaking signals are related to low ly- 
ing modes of the Dirac operator. The axial anomaly in- 
volves the topological charge of the gauge configuration, 
which is proportional to the net number of exactly chi- 
ral (zero-)modes via the Atiyah-Singer index theorem 121. 
The chiral condensate is associated with the density of 
the Dirac operator's low lying (but non-zero) modes [ij. 
The non-vanishing quark condensate indicates breaking 
of both symmetries. 



IV. GAUGE CONFIGURATIONS 



Ginsparg- Wilson equation and therefore exhibits better 
chiral properties than the simpler Wilson Dirac opera- 
tor while being less expensive by an order of magnitude 
- in terms of computation time - in comparison to the 
chirally exact overlap operator. 

We calculated up to the lowest 256 eigenmodes of the 
Dirac operator D and up to lowest the 512 eigenmodes of 
the hermitian operator D^ using ARPACK which is an 
implementation of the Arnoldi method to calculate part 
of the spectrum of arbitrary matrices [251 • 

The quark propagator S is determined by inverting 
the Dirac operator for a given source. Instead of us- 
ing point sources we use Jacobi smeared sources [2^, |23| 
which are approximately of Gaussian shape. Their shape 
was adjusted to a width of about 0.27 fm [T3|. The low 
mode contribution ^luisi^fc) to the quark propagator, see 
(H)), has to be multiplied with the same sources as the 
full propagator S in order to achieve the correct reduced 
propagators S'rcd5(fc)- 



V. MESONS 

We restrict ourselves to the study of isovectors, in par- 
ticular the chiral partners: 



The vector mesons p (J 



PC _ 1 — 



1 ) with interpo- 



lating fields u{x)^id{x) and u{x)^i'yid{x) and ai 
{J^'^' = 1^^) with interpolating field u(a;)7i75d(a;); 
in a chirally symmetric world the vector and the 
axial vector interpolator get mixed via the isospin 
axial transformations. 

• The pseudoscalar n {J^'^ — ^) with interpolat- 
ing fields u{x)^zd(x) and u{x)^i^zd{x). We also 
study the scalar ao {J^'~^ ~ 0++), u{x)d{x) which 
would get mixed with u{x)^Qd{x) via the \]{\)a 
transformation. 

(In the interpolators 74 denotes the Dirac matrix in Eu- 
clidean time direction.) 

We compute from the quark propagators meson propa- 
gators, projected to vanishing momentum and determine 
the hadron masses from a range of Euclidean time val- 
ues where the correlation function exhibits exponential 
decay. The final errors are statistical only and obtained 
by standard jackknife elimination sampling. 



For our analysis we used 161 gauge field configurations 
[ITI I18| of lattice size 16^ x 32; with the lattice spac- 
ing a = 0.144(1) fm this corresponds to a spatial size 
of 2.3 fm. The simulation includes with two degenerate 
flavors of light fermions and a corresponding pion mass 
of TOtt — 322(5) MeV. For the dynamical quarks of the 
configurations as well as for the valence quarks the so- 
called Chirally Improved Dirac operator [23, [2J] has been 
used. This operator is an approximate solution to the 



VI. RESULTS 



Low mode sector 



Figure [T] shows the integral over the distribution H{\n\) 
of the (real) eigenvalues of D5. The scale is set by the 
lattice spacing. There is a transition region up to roughly 
twice the size of the quark mass (for this simulation 
the unrenormalized mass calculated from the axial Ward 
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FIG. 1. The integrated eigenvalue density for the lowest 512 
(absolute) eigenvalues of D5. The eigenvalues are scaled ac- 
cording to the lattice spacing. The number on the upper axis 
indicate the values of fi where there are 16, 32, 64, 128, 256 
and 512 eigenvalues below that value. 




FIG. 2. Low mode contribution to the correlators for the 
J^ = 0~^ sector in comparison to the correlators from fuU 
propagators with interpolators (a) U'jsd and (b) u^4"/5d. The 
number of included modes is shown in the legend. 



identity is 15 MeV [l^) corresponding to 0(16) eigen- 
modes, as also observed in, e.g., |19l422l |. As will be seen 
below, this is in accordance with the behavior observed 
for the meson propagators. 

For the overlap operator the real eigenvalues corre- 
spond to exact chiral modes, the zero modes. This is no 
longer true for Wilson- type operators. There one may as- 
sociate zero modes with real eigenvalues, although there 
chirality is not unity. For the hermitian Dirac operator 
there is no simple method to identify these would-be zero 
modes and thus all we can say is that the lowest 8 modes 
include a significant number (if not all) of the would-be 
zero modes (instantons) . 

Before we construct meson correlators out of reduced 
quark propagators, let us first consider meson correlators 
approximated by the lowest k modes only, using propa- 
gators S'i,„5(fc), see (dl. 

In Fig. [5] we compare the pseudoscalar correlator us- 
ing standard full propagators to the correlators using 
only the lowest modes of the hermitian Dirac operator 
-D5. For the two pseudoscalar operators the exponential 
pion decay behavior sets it much earlier (at lower num- 
bers of eigenmodes) for the interpolator uj^d than for 
the other interpolator 117475^. Clearly the first one is 
stronger dominated by the low lying modes than the sec- 
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FIG. 3. Correlation functions for the reduced interpolators 
as compared to the correlators from full propagators. Top; 
J^ = 0""^ with interpolators (a) u'ysd, (b) 117475^. Middle: 
J^ = 1~~ with (c) u-jid, (d) U747id. Bottom: Reduced (e) 
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0++ (ud) and (f) J^ 



(u'yi'ysd). 



ond. The large time region is being well described by the 
low modes whereas the short time region - where excited 
states dominate - gets much slower saturated. Compar- 
ing with the result for an equivalent approximation for 
the non-hermitian Dirac operator (not shown here) we 
find that less eigenmodes of D5 are needed to obtain a 
similar quality of approximation of the correlators with 
full propagators. These results agree with the observa- 
tions in [aSlil- 



B. Removing the low mode sector 

Figure [3] shows the meson propagators for various 
stages of low mode removal, always in comparison with 
the full propagator and Fig. [S] combines the correspond- 
ing mass fits to the regions of exponential behavior. 

All mass values (except for the p) exhibit a strong de- 
pendence on the truncation of the lowest eigenmodes; 
from truncations levels of ^ 16 modes upwards (corre- 
sponding to quark masses of approximately 30 MeV) all 
mass values then follow a roughly parallel, rising behav- 
ior. The range of exponential behavior of the correlators 
shrinks, as can be seen in the log-plots in Fig. [31 

The effective mass plots (the local two-point approxi- 
mation of the derivative of the logarithm of the correla- 
tors) in Fig. |4] indicate the regions, where an exponential 
fit to the correlators has been done. We find that the 
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FIG. 4. Effective mass plots for the reduced interpolators as 
compared to the full propagators. For the notation (a-f) see 
Fig. El 
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FIG. 5. The masses of all considered mesons as a function of 
the reduced spectrum, subtracting the 0-512 lowest modes of 



fluctuation typically decreases with increasing reduction. 
This may be related to the relative importance of the 
noisy low lying modes in the quark propagators. 

In [2^ the parity-chiral group and the effect of symme- 
try breaking on the meson spectrum is discussed. E.g., 
whereas the U(l)^ breaking lifts the degeneracy between 
pion and ao (and between 77 and /o) the breaking of the 
chiral SU(2)i x 811(2)^? symmetry is related to the mass 
differences of pion and /o (and ao and rj). From Fig. 
[S] we find drastic sensitivity on low modes for both, the 



pion interpolator masses and the ao-mass. At low trun- 
cation levels the ao-mass rapidly drops; it does not drop 
down to the pion mass value. This might indicate some 
remnant of the anomaly breaking for the J — states. 

The pion interpolators exhibit a puzzling behavior. 
The classical pion interpolator uj^d quickly loses its ex- 
ponential behavior at larger (Euclidean) distances; only 
a more massive decay signal is observed at smaller dis- 
tances (Fig. 131). From trtmcation level 16 onwards we 
therefore do not exhibit mass values in Fig. [S] for that 
interpolator. A fit to the very small time slices gives a 
mass approaching the mass value from the second inter- 
polator U7475d with the pion quantum numbers, which 
couples due to PC AC (proportional to the quark mass). 

For the J^*-^ = 1^^ vector meson p there are two chi- 
ral representations, which correspond to the vector in- 
terpolator ujid and (Dirac-) tensor interpolator uj^'jid 
Their chiral partners [28] are the oi and the hi mesons, 
respectively. We did not determine the hi mass, since 
its interpolator includes disconnected graphs (it is an 
/ = state). There is no noticeable splitting between 
the two p-interpolators for all stages of truncation. We 
do find, however intriguing behavior comparing the p- 
mass with the ai result. Starting out quite differently 
for the full quark propagator the masses approach each 
other and are compatible with each other from trunca- 
tion level 8 onwards. This indicates restoration of the 
SU(2)l X SU(2)ii: symmetry for J = 1 states. The very 
fact that all three interpolators (vector, tensor and axial 
vector) give the same mass hints to the restoration of the 
SU(2)l X SU(2)fl^ X U(1)a symmetry for J = 1 states. 
The latter could be reliably concluded, however, only af- 
ter studying of the hi meson. 



VII. CONCLUSIONS 

The low lying eigenvalues of the Dirac operator are 
usually associated with chiral symmetry breaking. We 
have computed hadron propagators while removing in- 
creasingly more of the low lying eigenmodes of the Dirac 
operator. This allows us to study their influence on cer- 
tain hadron masses. Due to the relationship of the low 
eigensector with chiral symmetry breaking, this amounts 
to partially restoring chiral symmetry (in the valence 
quarks). 

We find drastic behavior for some meson interpolators 
when starting to remove low eigenmodes. At truncation 
level 16 the behavior saturates and then the mass values 
rise uniformly with roughly parallel slopes. The confine- 
ment properties remain intact, i.e., we still observe clear 
bound states for most of the studied isovector (scalar, ax- 
ial vector and vector) mesons. An exception is the pion, 
where no clear exponential decay of the correlation func- 
tion is seen in the uj^d interpolator, but a massive state 
is seen in the u'fi'y^d interpolator. The mass values of 
the vector meson chiral partners ai and p approach each 
other rapidly when 8 or more low modes are removed. 



We conclude that essential confinement properties re- 
main intact, even when the low eigenmodes of the Dirac 
operator are removed in the valence sector. Restoration 
of chiral symmetry is observed in that approximation. 
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